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The auxiliary field method is a powerful technique to obtain approximate closed-form energy for- 
mulas for eigenequations in quantum mechanics. Very good results can be obtained for Schrodinger 
and semirelativistic Hamiltonians with various potentials, even in the case of many-body problems. 
This method can also provide approximate eigenstates in terms of well known wavefunctions, for 
instance harmonic oscillator or hydrogen-like states, but with a characteristic size which depends 
on quantum numbers. In this paper, we consider two-body Schrodinger equations with linear, log- 
arithmic and exponential potentials and show that analytical approximations of the corresponding 
eigenstates can be obtained with the auxiliary field method, with a very good accuracy in some 
, cases. 

q ■ PACS numbers: 03.65.Ge 

(N : 
c : 

1—5 1 I. INTRODUCTION 

m : 

A considerable amount of works have been already devoted to the computation of closed-form solutions of eigenequa- 
tions in quantum mechanics, especially in bound states problems. Apart from its intrinsic mathematical interest, 
Oh finding analytical formulas is always useful in physics, for example to obtain informations about the dependence of 
+L> . observables on the parameters of the Hamiltonian and on the quantum numbers. 

Many methods exist in the literature to find approximate analytical solutions of eigenequations for bound state 
problems: perturbation theory, Wentzel-Kramers-Brillouin (WKB) method, large- N expansion, variational method, 
etc. (see for instance Refs. [2, HJ ). Recently, we have proposed a new technique, that we called the auxiliary field 
i—i ] method (AFM), to compute bound states of a Schrodinger equation Q. This method is based on the introduction 
of auxiliary fields - also known as einbein fields [3] _ into the original Hamiltonian and can lead to approximate 
closed- form results for various potentials [B-0|- It has been extended to treat relativistic Hamiltonians Q and even 
many-body problem s M . It is worth mentioning that we have shown in a previous work [Tj| that the AFM and the 
envelope theory [lll - ll6| are, to a large extent, two equivalent approaches. For instance, pure powers potential have 



been treated in detail in Refs. and analytical results have been obtained for the Coulomb- plus- linear interaction 

7-H ' in Refs. dEl. 

Since the principle of the AFM is to replace the studied Hamiltonian H by a Hamiltonian H for which analytical 
solutions are known, this method can also provide approximations for the genuine eigenfunctions of H in terms of the 
eigenfunctions of H. The trial solutions are then characterized by a size which depends on the quantum numbers. 
^ — 1 ■ The purpose of this paper is to test the quality of the trial wavefunctions given by the AFM by computing several 
observables, like mean values of powers of the operators position and square momentum. We consider two-body 
• *h , Schrodinger equations with three different potentials, relevant in particular for hadronic physics: linear, logarithmic, 
and exponential ones. We focus particularly on the case of the linear potential because analytical solutions are known 
^ for S-eigenstates. 

Our paper is organized as follows. The general principles of the AFM are recalled in Sec. [II] The method is applied 
in Sec. IHII to the case of the linear potential using two different approaches. Logarithmic and exponential potentials 
are studied in Sees. HVl and IVl respectively. Finally we sum up our results in Sec. IVII In order to be complete, a lot 
of mathematical expressions are given in this paper but, to avoid that technical details hamper the main text, they 
are relegated in appendices. Formulas to compute some observables are given in App. |2] Properties of Airy function, 
hydrogen-like system and harmonic oscillator are respectively reminded in App. [B] [C] and [D] The overlap of some 
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dilated eigenstates are detailed in App. |EJ Finally, the inverse of a special function is also given in App. [F] 

II. THE AUXILIARY FIELD METHOD 

In order for the present paper to be self-contained, we recall here the main points of the AFM and refer the 
reader to Refs. 0, for more details. The aim of this method is to find an analytical approximate solution of the 
eigenequation H\ip) — E\ip), with 

H = T(p 2 ) + V(r), (1) 
where V(r) is a central potential. Let us consider a new potential 

V(r,y) = yP{r) + V (I{y)) - »P (/(»)) (2) 

for the Hamiltonian 

H(v)=T(p 2 ) + V(r,i>), (3) 

which depends on the auxiliary field v. Up to a dimensioned factor, P(x) plays the role of another potential and I(x) 
is the inverse of a function K(x); explicitly, these functions are defined by 

I{x)=K-\x) and K(x) = ^-, (4) 

the prime denoting the derivative with respect to x. A proper elimination of this auxiliary field v by the following 
constraint SoH(i))\~ = ~ o = leads to the solution Vq = K(r). The original Hamiltonian ((T|) can then be recovered 

since V(r,0 ) = V(r). 

The idea of the AFM is then to replace this operator v by a real parameter v. If the potential P(r) is well chosen 
(for instance P(r) = —l/r or r 2 ), the eigenequation 

H A \1>(v)) =E A {v)\^{v)) (5) 

with Ha = T(p 2 ) + vP(r) is analytically solvable for a nonrelativistic kinetic energy operator, as well as the one 
associated with the Hamiltonian 

H{u)=T{p 2 )+V{r,v). (6) 

Let us remark that V(r, v) is of the form vP(r) + C(y) where v and C(y) are constants. Using (JSJ) , an eigenvalue 
E(v) of H{y) is given by 

E{v)=E A {v) + V{I{v))-vP{I{v)). (7) 

The approximate eigenvalues and eigenstates are eventually given by E{uq) and \rft(yo)) respectively, with vq such 
that E{vq) is extremal, i.e. 

WU = o- (8) 

The value of vq depends on the quantum numbers of the state considered and should be denoted more precisely as 
fo(n, I). It is worth noting that if P(r) = V(r), the method gives the exact results. Despite the fact that the auxiliary 
field is replaced by a real parameter, we keep the name "auxiliary field method" used in previous studies [l7| . In these 
papers, the technique was applied to get rid of the square root operator appearing in semirelativistic Hamiltonians 
(see also Refs. (Hi). 

Once Vq is known, the constant C(vq) can be computed and is such that V(tq, vq) = V{tq) and V'(ro, vq) — V'(tq), 
where tq = I{vq). It can be shown that Q 

(^)|P(r)|^))=P(r ). (9) 

This means that tq is a kind of "average point" for the potential P(r). Using this last relation with the definitions 
above, we get 

(V>(z/ )|J(*>oM^o)> = JOo) with J(x) = P(I(x)). (10) 
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So, our method can actually be considered as a "mean field approximation" with respect to a particular auxiliary held 
which is introduced to simplify the calculations: v§ is the mean value of the operator Oq — K{r) through a function 
J which can be quite simple. For example, J(x) = x if V(x) = P(x) 2 /2 + Vb where Vq is a constant. 

Since the potential V(r, v$) is tangent to the potential V(r) at r — ro, the approximation E(i/q) is an upper (lower) 
bound of the exact energy if V(r, Vq) > V(r) (V(r, vo) < V(r)) for all values of r [l(J. Equivalently, a function (7(2) can 
be defined by V(x) = g(P(x)). It can then be shown that, if g(x) is a concave (convex) function, the approximation 
E(vq) is an upper (lower) bound of the exact energy. This property has been demonstrated in the framework of the 
envelope theory [llTllq ]. but can be applied as well to the AFM [l(|. Several examples will be presented below. The 
knowledge of a lower and an upper bounds for an eigenstate is a first technique to estimate the accuracy of the AFM. 
It has also been shown that Q 

E{ VQ ) ~ = v(r ) - (VMinOIV-M). (ii) 

The right-hand side of this equation is the difference between the value of potential V computed at the average point 
ro and the average of this potential for the AFM state \ip(vo)) considered here as trial state. In some favorable cases 
(the trial state is a ground state for instance), E < {ip(i/o)\H\ip{yo)) and a bound on the error can be computed by 

E(uo) -E> V(r ) - Mn))|V(r)Mi«,)). (12) 

If the mean value of V can be computed analytically, this constitutes a second procedure to estimate the accuracy 
of the AFM. Several examples of this calculation are presented in Ref. [3j. At last, the eigenstates of a Hamiltonian 
of type ((T|) can be solved numerically with an arbitrary precision. So, as a third possibility, a direct comparison 
with the AFM results is always possible. So, one can wonder why to use the AFM? The interest of this method is 
mainly to obtain analytical information about the whole spectra (dependence of eigenenergies on the parameters of 
the Hamiltonian and on the quantum numbers) , without necessarily searching a very high accuracy. 

The AFM is completely general and a priori valid for any potential V(r). Nevertheless, in order to get closed-form 
expressions, we must fulfill two conditions: first, to be able to determine the function I(x) defined by §4§ and, second, 
to be able to find vq defined by ([5]) and to calculate the corresponding value E(v$) in an analytical way. Our previous 
works showed this can be managed for a lot of interesting physical problems and that very accurate approximate 
values can be obtained for the energy. 

Moreover, the state which is an eigenstate of H(vq), is an approximation of a genuine eigenstate of H. 

The shape of the corresponding wavefunction (r\^{yo)) depends on the quantum numbers via the parameter vq. In 
the following sections, the quality of this approximation will be tested for the Schrodinger equation with different 
potentials by using the AFM with two different functions P{r) — r 2 and — 1/r. 

As it is shown in Ref. Q, the same energy formula E(N) is obtained for the choice P(r) = sgn(A)r A with A > —2, 
but with different forms for the global quantum number N. N = 2n + 1 + 3/2 for A = 2, and N = n + 1 + 1 for A = — 1. 
For other values of A, The form of N is not exactly analytically known. 



III. THE LINEAR POTENTIAL 



Let us denote m the mass of the particle for a one-body problem or the reduced mass for a two-body system. 
When T(p 2 ) = p 2 /(2m) and V(r) = ar with a positive, the Hamiltonian H can be solved analytically for a vanishing 
ang ular momentum I (see App. IB}, When I ^ 0, numerical results have been obtained with two different methods 
[l8l li~9l|. In order to simplify the notation, we will denote \n,l) an eigenstate of H and \n) = \n, 0), with (r\n) given 
by (IB3|) . The corresponding energies will be denoted E n< i — (n, l\H\n, I) and E n = E n ,o — (n\H\n). As we need a 
Hamiltonian H with a central potential which is completely solvable to apply the AFM, we can only use in practice 
a hydrogen-like system (P(r) = —1/r) or a harmonic oscillator (P(r) = r 2 ). These two Hamiltonians will be used 
respectively in Sec. IIII Al and IIIIBI As a linear potential seems closer to r 2 than —1/r, we can expect that the use 
of a harmonic oscillator to start the AFM will give better results. The scaling properties of the approximate AFM 
solutions and the exact solutions are the same [3, Q, so we can set m = 1/2 and a = 1 to lighten the notations and 
to match the results obtained in Ref. [3(. So, the following Hamiltonian is considered 

H = p 2 + r. (13) 

We denote e n ,i the approximate AFM energies which are given by Q 

^ = 4^ 2/3 - (14) 
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Exact energies, given by (|B2[) . reduces to E n — —a n , where a n is the (n+ l) th zero of the Airy functions Ai. A simple 
approximation of E n can be obtained using the expansion (|B4I) at the first order 



x 2/3 / n\ 2/3 / oX 2/3 



E n = -a n n ^—j ~ 2,811 l n+ lj ' (15) 

It is worth noting that the linear potential is not only a toy model to test the AFM method. Effective theories of the 
quantum chromodynamics have proved that it is agood interaction to take into account the confinement of quarks 
or gluons in potential models of hadronic physics [13, H3| ■ 



A. AFM with P(r) = -1/r 

One can ask whether it is possible to obtain good approximations for the solutions of a Schrodinger equation with 
a linear potential by means of hydrogen-like eigenfunctions. This will be examined in this section. Using the results 
of Ref. j3j for P(r) = —1/r, we find r = y/vo with v a = 2 2 / 3 (n + / + l) 4 / 3 . Exact eigenstates are then approximated 
by AFM eigenstates which are, in this case, hydrogen-like states (|C3[) with 

v (n + l + l)^ 3 

r > = Y = W 3 ■ (16) 

Such states are denoted |Hy;n, I) and |Hy;n) = |Hy;n, 0). Using (|C6[) and results above, it can be shown that © 
is satisfied, with P(r ) = —1/tq = — 2~ 1//3 (n + I + 1)~ 2 / 3 . It is also worth noting that (fl"0)) gives (l/Vt^ = 1/y/Uo. 
We denote and e^y = the approximated energies which are given by (fT4| with N = n + I + 1 3J. Since 
V(r, i^o) — V(r) — — (r — tq) 2 /r < 0, are lower bounds of the exact energies. This can also be determined with the 
function g defined in Sec. HU In this case, g(y) = —1/y with y < 0. The function g"(y) — —2/y 3 being positive, g is 
convex as expected. 

The quantum number dependence of the scaling parameter 77 corrects partly the difference between the shapes of 
(r\n, I) and (r|Hy; n, I). Consequently, (Hy; n, l\Hy; n, I') = 8w because of the orthogonality of the spherical harmonics, 
but (Hy; n, Z|Hy; n', I) 7^ 8 nn i. Using the definition (jEll) . we find 



'n' + 1 + 1 



4/3N 



(Hy; n, Z|Hy; n> , I) = F«*. , « ^ ^^y^j j (17) 

with F Hy given by (|E3|) . TableUgives some values of |(Hy; n, Z|Hy; n', I) | 2 . We can see that the overlap is not negligible 
for n close to n', but it decreases rapidly with \n — n'\. The situation improves when I increases: |(Hy; 0, ^|Hy; 1, l)\ 2 = 
0.43, 0.29, 0.22, 0.18, 0.15, 0.13 for 1 = -> 5. 



TABLE I: Results for P n ,„/,! 


= |(Hy;n,Z|Hy;n',Z>| 2 . 


Values for / = (I = 1) are 


given in the lower-left (upper-ri; 


5M) triangle 


of the Table. P n>n /,j = P n >,n,t 


and P„, n ,i = 1. 










n = 


1 


2 


3 


n' = 


1 


0.29 


0.028 


0.0039 


1 


0.43 


1 


0.36 


0.036 


2 


0.055 


0.43 


1 


0.39 


3 


0.0097 


0.049 


0.43 


1 



1. Results for 1 = 

It is interesting to compare e^ y with (TTS"]) 

^ y -i(- + 1 ) 2/3 - 1 - 89 °("+ 1 ) 2/3 - ("J 

The ratio e^ y /£" n is respectively equal to 0.808, 0.734, 0.712, for n = 0, 1, 2 and tends rapidly toward the asymptotic 
value S 1 / 3 /^ 2 ' 3 ~ 0.672. As expected, these ratios are smaller than 1 since e^ y are lower bounds. Two wavefunctions 
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are given in Fig. [TJ We can see that the differences between exact (r\n) and AFM (r|Hy;n) wavefunctions can be 
large. The overlap |(n|Hy;n)| 2 between these wavefunctions can be computed numerically with a high accuracy. 
We find respectively the values 0.934, 0.664, 0.298 fo r n = 0,1,2, showing a rapid decreasing of the overlap. It is 
worth noting that, asymptotically, (r\n) oc exp (— | \/2mar 3 / 2 ) while (r|Hy;n) is characterized by an exponential 
tail. Nevertheless, if an observable is not too sensitive to the large r behavior, this discrepancy will not spoil its mean 
value. 




FIG. 1: Normalized wavefunctions (r\n) (solid line) and (r|Hy;n) (dashed line) for n — (left) and 1 (right). 



The various observables \ip n fi(0)\ 2 , (r k ) and (p k ) computed with the AFM states can be obtained using formu- 
las (|C4[) . (|C6[) and (|C7[) for I = with the parameter r] given by (HI))) . Results are summed up in Table HB A direct 
comparison between the structure of exact and AFM observables can be obtained if we remind that the exact ones 
depend on \a n \ which can be well approximated by j3 n (see (|B4|) ). Let us look in detail only at the mean value (r). 
For the exact and AFM solutions, we have respectively 

2|on| /2tt 2 \ 1/3 / 3\ 2/3 / 3^ 2/3 



(n\r\n) = « [~ ) { n + 4 ) «l-874ln+-l , (19) 

(Hy;n|r|Hy;n) = JL(n + l) 2 / 3 w 1.890 (n + 1) 2/3 . (20) 

Some observables like |-0 nj o(O)| 2 and (p 4 ) are very badly reproduced. Others can be obtained with a quite reasonable 
accuracy. Despite the fact that exact and AFM wavefunctions differ strongly when n increases, their sizes stay similar. 

Since (Hy; n|p 2 |Hy; n) and (Hy; n|r|Hy; n) are known, it is possible to compute analytically e^ y = (Hy; n|iT|Hy; n). 
|Hy; n), which is an exact eigenstate of the approximate Hamiltonian H , is then considered as a trial eigenstate for the 
genuine Hamiltonian H. The results, given in Table [TTl show that e^ y is a better approximation than e^ y . Because 
of the Ritz theorem, e^ y > E , but the variational character of other values e^ y cannot be guaranteed, in contrast to 
the values e^ y . 



TABLE II: Various observables Q computed with the AFM and P(r) = — 1/r. The ratios e^ y /E n given in the last row can be 
compared with the ratios eJJ y /_E n given in the penultimate row. 





(Hy;n|Q|Hy;n) 






(Hy;n|Q|Hy;n)/(n|Q 


n) 


(Q) 


(2m = o = l) 


n = 


n = 1 


n = 2 


n > 1 


l^,o(o) | a 


n+'l 
2tt 


2 


4 


6 


2(n + l) 


(r) 


3 ( n+ l)2/3 
2 2/3 


1.212 


1.101 


1.068 


1.009 + 2^ + 0(4,) 


(r 2 ) 


5n +10n + 6 
2 1 /3( n + l)2/3 


1.633 


1.178 


1.080 


0.942 + M= + (J T ) 


(r 3 ) 


5(7n 2 + 14n + 12) 
4 


2.392 


1.303 


1.099 


0.862 + + G (J T ) 


(r 4 ) 


3(21n 4 +84n 3 + 161n 2 + 154n + 60) 
25/3( n + l)4/3 

(n + 1) 2 / 3 
2 2/3 


3.752 


1.517 


1.148 


0.782 + 2^22 +0(4,) 


(P 2 ) 


0.808 


0.734 


0.712 


0.672 + 2=12. + 0(-*j) 


</> 


(8n + 5)(n + l) 4 /3 
2 4/3 


1.815 


3.890 


5.916 


2.009 n+ 1.926 + 0(i) 


= ^ y 


(P 2 ) + (r) 


1.078 


0.978 


0.949 


0.896 + 2i= + o(^) 






0.808 


0.734 


0.712 


0.672 + 2i= + o(^) 
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2. General results 

The behavior of observables computed with the AFM is similar for values of I = 0,1,2. We do not expect strong 
deviations for larger values of I. This is illustrated with some typical results gathered in Table IPLT1 Observables are 
generally not very well reproduced, but it is expected since eigenstates for a linear potential are very different from 
eigenstates for a Coulomb potential. Actually, the agreement is not catastrophic, except for \ip n> o(0) | 2 and (p 4 ) as 
mentioned above. It is even surprising that the AFM with P(r) = — 1/r could give energies and some observables for 
a linear potential with a quite reasonable accuracy. Moreover, lower bounds of the energies are obtained. This can 
be useful for some calculations (see Sec. IIIIC|) . 



TABLE III: Ratios between the AFM results (energies and (r)) with P(r) — —1/r and the exact results, for several quantum 
number sets (n,l). 



1 


n = 


n = 1 


n = 2 


n = 3 


n = 4 


n — 5 



















0.808 


0.734 


0.712 


0.702 


0.696 


0.692 


1 


0.893 


0.805 


0.767 


0.746 


0.733 


0.724 


2 


0.925 


0.846 


0.805 


0.779 


0.762 


0.749 


(Hy; n\r 


\Ry; n)/(n\r\n) 















1.212 


1.101 


1.068 


1.053 


1.043 


1.037 


1 


1.116 


1.118 


1.103 


1.089 


1.079 


1.071 


2 


1.080 


1.110 


1.110 


1.104 


1.096 


1.089 



B. AFM with P(r) = r 2 

Since the quadratic potential is closer to the linear potential than a Coulomb one, one can expect better results with 
harmonic oscillator wavefunctions. This will be examined in this section. Using the results of Ref. 3] for P(r) = r 2 , 
we find ro = l/(2vo) with vq — 2 _4 / 3 (2n + Z + 3/2)~ 2 / 3 . Exact eigenstates are then approximated by AFM eigenstates 
which are, in this case, harmonic oscillator states (|D3[) with 

A = v\l A = 2" 1/3 (2n + I + 3/2)" 1/6 . (21) 

Such states are denoted |HO;n, I) and |HO;n) = |HO;n, 0). Using (ID6[) and results above, it can be shown that © 
is satisfied, with P(r ) = r 2 = 2 2 / 3 (2n + 1 + 3/2) 4 / 3 . It is also worth noting that ([ID) gives (l/t> 2 ) = 1/Vq. This is in 
agreement with (20) in Ref. § (an auxiliary field <j>= l/v\s used in this last reference). We denote e^f and = 

the approximated energies which are given by (HD with N = 2n+l+3/2 Q. Since V(r, Vq )— V(r) = (r— r ) 2 /(2r ) > 0, 
are upper bounds of the exact energies. In this case, g(y) = y/y with y > 0. The function g"(y) = — l/(4y 3 / 2 ) 
being negative, g is concave as expected. 

The scaling parameter A depends on the quantum numbers. Consequently, (HO; n, Z|HO; n, I'} = Sw because of the 
orthogonality of the spherical harmonics, but (HO; n, l\HO; n', I) 7^ S nn /. Using the definition (|E1|) . we find 

(HO; n , i |HO; n ',0^«.,((^|i|) ,/6 ) (22) 

with F HO given by (|E4|) . Table ITVl gives some values of |(HO; n, /|HO; n', l)\ 2 . We can see that the overlap is always 
small and decreases rapidly with \n — n'\. The situation is even better when I increases: |(HO; 0, /|HO; 1, l)\ 2 — 0.029, 
0.023, 0.019, 0.017, 0.014, 0.013 for I = 5. 

1. Results for 1 = 

Let us look at the energies e; 



e»° = 3 



2/3 



(23) 
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TABLE IV: Results for P„ 


B /,, = |(HO;n,Z|HO;n',0| 2 . 


Values for I = (Z = 1) are 


given in the lower-left (upper-rif 


;ht) triangle 


of the Table. P nt „ij = P n i 


n> i and P„, n ,i = 1. 










n = 


1 


2 


3 


n' = 


1 


0.023 


0.0026 


0.00039 


1 


0.029 


1 


0.026 


0.0027 


2 


0.0036 


0.027 


1 


0.026 


3 


0.00064 


0.0031 


0.027 


1 



By comparing with (|15[) . we can see immediately that the situation is more favorable than in the previous case. The 
ratio e^°/E n is respectively equal to 1.059, 1.066, 1.067, for n = 0, 1,2. The asymptotic value 3 1 ^ 3 2 2 ^/t: 2 ^ m 1.067 
is rapidly approached. As expected, these ratios are greater than 1 since are upper bounds. Two wavefunctions 
are given in Fig. [51 We can see that the differences between exact (r\n) and AFM (r|HO;n) wavefunctions are quite 
small. The overlap |(n|HO; n)| 2 between these wavefunctions can be computed numerically with a high accuracy. 
We find respectively the values 0.997, 0.979, 0.951 for n — 0,1,2. A value below 0.75 is reached for n = 6 and 
below 0.25 for n = 14. A wavefunction (r|HO;n) is characterized by an Gaussian tail, while we have asymptotically 
(r|n) oc exp (— |%/2r7ia r 3 / 2 ) . Again, if an observable is not too sensitive to the large r behavior, this discrepancy will 
not spoil its mean value. 




FIG. 2: Normalized wavefunctions (r\n) (solid line) and (r|HO;n) (dashed line) for n — (left) and 1 (right). 



The various observables |VVt.o(0)| 2 , (r ) and (p k ) computed with the AFM states can be obtained using formu- 
las (|D4I) . (|D6[) and (|D7|) for I = with the parameter A given by (12T1) . Results are summed up in Table fVl Again, a 
direct comparison between the structure of exact and AFM observables can be obtained by using (3 n (see (|B4I) ) instead 
of \a n \. Let us look in detail only at the mean value (r 2 ). For the exact and AFM solutions, we have respectively 



. . 2| . 8|a„| 2 2 5 / 3 3 1 /3^/3 / gX 4/3 j ,4/3 

(n\r 2 \n) = « g \n+-j w 4.214 in + - J , (24) 



3 



4/3 



(HO;n|r 2 |HO;n) = 4 In + - J . (25) 

In contrast with the previous case, all observables are very well reproduced. This is also the case for n 3> 1, while the 
overlap |(n|HO;n)| 2 tends towards very small values in this limit. This is due to the fact that some observables are 
not very sensitive to the details of the wavefunctions and that the sizes of exact and AFM states stay similar. 

In this case, it is also possible to compute analytically = (HO; n\H |HO; n). The results, given in Table fVl show 
that is a better approximation than e^°. Because of the Ritz theorem, > Eq, but the variational character 
of other values cannot be guaranteed, in contrast to the values e^°. 



2. General results 



As in the previous case, the behavior of observables computed with the AFM is similar for values of I — 0,1,2. 
We do not expect strong deviations for larger values of I. Some typical results are presented in table IVT1 Agreement 
between AFM and exact results are very good, much better than for the previous case. This is expected, since 
eigenstates for a quadratic potential are closer to eigenstates for a linear potential than eigenstates for a Coulomb 
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TABLE V: Various observables Q computed with the AFM and P(r) — r 2 . The ratios e^°/E n given in the last row can be 
compared with the ratios /E n given in the penultimate row. 





(HO;7i|Q|HO;n) 




(ftO;n\Q\nO;n)/(n\Q\n) 




Q 


(2m = a = 1) 


n = 


n = 1 n — 2 


n > 1 



|^,o (0) | 2 ^nyan+e ^ 0^05 O902 0.900 + + O (4,) 



/5 - 



(r) 4(8n+b) 0.976 0.964 0.962 

(r 2 ) 4(n + 3/4) 4/3 0.935 0.946 0.948 0.949 - + Q 

(r 3 ) 8^(4n+3)»/»r(n+ 3 /2) 0.881 0.934 0.941 0.946 -^ + 0(^3 



r 4 } §(8n + 6) 2/3 (8n 2 + 12n + 5) 0.819 0.918 0.934 0.946 -^ + 0' — 

o 2 ) (n + 3/4) 2/3 1.059 1.066 1.067 1.067 - MSi + O ' 



7/3 

- 4 ) 3( f" 2+12 "t 5) 1.039 0.966 0.956 0.949 + 2^ + Q ( -k* 

1 4(8n + 6) 2 / 3 n* \ n 7 ' 3 j 

(H) = (p 2 ) + (r) 1.004 0.998 0.997 0.996 + ^ + q 

1.059 1.066 1.067 1.067 - + Q / l 



potential. The quantum number dependence of the scaling parameter A corrects, much better than 77, the difference 
between the shapes of AFM and exact eigenstates. 



TABLE VI: Ratios between the AFM results (energies and (r)) with P(r) = r 2 and the exact results, for several quantum 
number sets (n,l). 



1 


n = 


n = 1 


n = 2 


n = 3 


n = 4 


n = 5 





1.059 


1.066 


1.067 


1.067 


1.067 


1.067 


1 


1.036 


1.055 


1.060 


1.063 


1.064 


1.065 


2 


1.026 


1.046 


1.054 


1.058 


1.061 


1.062 


(HO;n|r|HO; 


n) / (n\r\n) 















0.976 


0.964 


0.962 


0.962 


0.961 


0.961 


1 


0.985 


0.972 


0.968 


0.965 


0.964 


0.963 


2 


0.990 


0.978 


0.972 


0.969 


0.967 


0.966 



C. Supplementary topics 

For the linear potential, the energy formula obtained with the AFM is e(N) = 3(A r /2) 2 / 3 , with N = n + I + 1 for 
P(r) = — 1/r and with N — 2n + I + 3/2 for P(r) = r 2 . It has been shown that the accuracy of this energy formula 
can be greatly improved by changing the structure of N Q. A very convenient form for the linear potential is (see 
(64) and (69) in Ref. 0) 

N = JL n + l + ^L i (26) 

which gives 

This formula can be compared with (|T5|) . The scaling parameters of the AFM eigenstates (see (TTBl and ([2~Tj)L as well 
as the observables, depend explicitly on this global quantum number N. So, one can ask whether it is also possible to 
improve the accuracy of these observables by using the formula (|26l) instead of the natural values of N. Unfortunately, 
this replacement spoils the values of some observables. We also tried to use a general form N — An + B I + C with 
A, B, C as free parameters. In so doing, some observables can be improved, but, in the same time, others are 
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dramatically spoilt. No unique set of values was found to improve globally the observables studied. Even in this case, 
no guarantee could be obtained about the behavior of others observables, and no systematic technique could be used 
to find the optimum parameters. So, if an improvement of the AFM energy can generally be obtained 0, 0-0], n o 
such improvement seems possible for a set of observables. 

The Eckart bound Be (see App.[A} is a lower bound on the overlap |(</>|0)| 2 between the exact ground state |0) and 
a trial state \<j>). In our case, it needs the computation of (<fi\H\(j>), with \<j>) — |Hy;0) or |HO;0), and the knowledge 
of the exact energies for the ground and the first excited states. But all these values are known (see above). From a 
practical point of view, it is generally only possible to compute the less accurate bound B' E defined in App. IA1 since 
only lower and upper bounds on energies must be known. Such quantities can be computed with the AFM (e ^ are 
lower bounds and are upper bounds), and the bound B' E can also be determined. Results, given in Table [yTTT 
show that the bound B' E , which is the easiest computable, cannot give relevant information in this case. On the 
contrary, the knowledge of the exact ground and first excited energies allows a good estimation of the overlap |(0|O)| 2 , 
without the necessity to compute the exact ground state |0). 



TABLE VII: Overlap |{<^|0)| 2 , computed numerically, between a AFM state \<f>) and the exact ground state |0) compared with 
the lower bounds Be and B' E (see App. [A}. 



w 


l(0|O)| 2 


B E 


B' E 


|Hy;0> 


0.934 


0.896 


0.195 


|HO;0) 


0.997 


0.995 


0.265 



IV. THE LOGARITHMIC POTENTIAL 



The logarith mic p otential, which has the particularity that its relative spectrum is independent of the mass of 
the particle [13, l2l| . can be used to simulate the confinement for heavy mesons. As the scaling properties of the 
approximate AFM solutions and the exact solutions are the same 0,0], we can work with reduced variables. In order 
to match the results obtained in Ref. Q, the following Hamiltonian is considered 

p 2 

H=^- + lnr (28) 

Accurate numerical solutions have been obtained with two different methods [3 [l^] . The notations of the previous 
section are reused to denote energies and eigenstates. Using the virial theorem, it can be shown that (n,l\p 2 \n, I) = 2 
for any eigenstate. The logarithmic potential can be considered as the limiting case of a power potential when 
the exponent tends toward zero pj. [20j. So, in some sense, this interaction is "midway" between parabolic and 
Coulomb potentials. It is then interesting to consider again both approximations given by P(r) = — 1/r and r 2 . The 
approximated AFM energies are given by Q 

e»,J - ^ (Jl N ) ■ ( 29 ) 



A. AFM with P(r) = -1/r 

Using the results of Ref. [|| for P(r) = we find ro = v$ with v$ — N j\[2 and N = n+ 1 + 1. AFM approximate 

eigenstates are hydrogen-like states (IC3|) with 

T) = 2iv Q = V2N. (30) 

In this case, g(y) = — ln(— y) with y < 0. The function g"{y) = 1/y 2 being positive, g is convex and e n ^ are lower 
bounds of the exact energies. 

Mean values (Hy; n, Z |r 2 |Hy; n, I) and (Hy; n, l\p 2 \Ry; n, I) are computed using (|C6[) and (|C7[) with ([SU]). In particular 
we have (Hy; n, l\p 2 \Hy; n, I) — 2. The AFM result is the exact one in this case. 
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B. AFM with P(r) = r 2 

Using the results of Ref. [3] for P(r) = r 2 , we find r = 1/v^o with ^ = V^ 2 and TV = 2n + I + 3/2. AFM 
approximate eigenstates are harmonic oscillator states (|D3I) with 

A=(4^) 1/4 = y|. (31) 

In this case, g(y) = \n(y/y) with y > 0. The function = — l/(2y 2 ) being negative, g is concave and e^y* are 

upper bounds of the exact energies. 

Mean values (HO; n, Z|r 2 |HO; n, I) and (HO; n, l\p 2 \RO; n, I) are computed using (fD6|) and (|D7| with f3lj). In 
particular we have (HO; n, Z |p 2 |HO; n, I) — 2. The AFM result is also the exact one in this case. 



C. Results 



Results concerning some observables derived from Hamiltonian (j2"5)l are gathered in Table IVIIII One can see that 
e y , and enclose the exact energy E n< i, as expected. Lower and upper bounds are not very good, mainly for the 
ground state. They approach the exact values when n or/and I quantum numbers increase, and the relative error 
decreases below 10% for the excited levels presented. The exact result is very close to the average between the lower 
and the upper bounds. It is worth noting that values of energies can be greatly improved (see Ref. @) but the 
variational character cannot be guaranteed anymore. 

As mentioned, AFM reproduces the exact value of {p 2 } for both HO and Hy wave functions. For unknown reasons, 
AFM has sometimes very astonishing virtues like this one. The AFM (r 2 ) values do not exhibit particular trends as 
function of the quantum numbers. They fluctuate around exact values but are maintained within satisfactory limits 
(less than 10%) for both HO and Hy wave functions. 

The overlap between AFM and exact wave functions is rather good for n = but deteriorates rapidly with increasing 
values of n, the role of / being much less important (a factor very important for the overlap is the number of nodes 
in the wave function). Here again HO and Hy wave functions give very similar results and behaviors. 

It may seem strange that the quality of the energies increases sensitively while, in the same time, the quality of 
the overlap decreases drastically. The same conclusion is valid for other observables. This is due to the fact that the 
overlap is very dependent on the details of the wave function, but that the considered observables depends essentially 
on the size of the wave function which is very similar for the exact and AFM wave functions. Contrary to the linear 
potential, the quality of the AFM results is quite similar whatever the choice of P(r) indicating that the logarithmic 
potential, being midway the parabolic and the Coulomb potentials, can be approached equivalcntly with HO or Hy 
type wave functions. We have studied also other observables, as for the linear case, but the corresponding results 
leading essentially to the same conclusions, they are not presented here. 



V. THE EXPONENTIAL POTENTIAL 



The exponential potential is very different from the parabolic, linear, and logarithmic ones. It is not a confining 
interaction since continuum states are allowed. So, it is closer to the Coulomb interaction, with the difference that 
only a finite number of bound states exists. Let us note that the exponential potential can be used to simulate 
screening color effects in mesons [22j. In order to match the results obtained in Ref. @, the following Hamiltonian 
written in reduced variables is considered 

H = p 2 -ke- r . (32) 

Accurate numerical solutions have been obtained with two different methods [HI, [l| . All bound states allowed have 
a negative energy. When the last excited states is characterized by an energy very close to zero, the corresponding 
wavefunction can have a very large extension. The notations of the previous section are reused to denote energies and 
eigenstates. Again, we consider both approximations given by P[r) = — 1/r and r 2 . The approximated energies are 
given by 



-ke 3W °^ (l+^W (T)\ 



= -A' T ' I 7 — + %-^\ with T = — - ( —) ' , (33) 

\Wq(T) 3 2Wq(T) 2 J 3 V k J V 1 
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TABLE VIII: Some AFM results for the Hamiltonian 
{X;n,l\p 2 \X;n,l) / (n,l\p 2 \n,l), and overlap |(n, l\X; n, 1} \ 2 . 


R(E) = eljE nM R(r 2 )-- 


= (X;n,l\r 2 \X;n, l)/(n 


,,l\r 2 \n,l),R{p 2 ) = 


For each set (l,n), results 


are given for P(r) = 


r 2 (X=HO) in the 


first line and for P(r) = 


— 1/r (X=Hy) in the second line. 








I n 


X 


R(E) 


R(r 2 ) 


R(p 2 ) 


\(n,l\X;n,l)\ 2 





HO 


1.591 


0.938 


1 


0.989 




Hv 




1.251 


1 


0.983 


1 


HO 


1.218 


1.051 


1 


0.766 




Hv 


U. ( oo 


0.900 


1 


0.771 


2 


HO 


1.164 


1.074 


1 


0.447 




Hv 


fl 784 


0.817 


1 


0.408 


1 


HO 


1.128 


0.959 


1 


0.987 




Hv 


u.oyo 


1.151 


1 


0.984 


1 


HO 


1.137 


1.014 


1 


0.792 




Hy 


0.859 


1.002 


1 


0.835 


2 


HO 


1.127 


1.041 


1 


0.470 




Hy 


0.856 


0.924 


1 


0.522 


2 


HO 


1.065 


0.969 


1 


0.987 




Hy 


0.948 


1.107 


1 


0.984 


1 


HO 


1.097 


0.998 


1 


0.815 




Hy 


0.909 


1.039 


1 


0.866 


2 


HO 


1.101 


1.020 


1 


0.499 




Hy 


0.895 


0.980 


1 


0.594 



where Wo is one branch of the Lambert function 6]. In this last paper, the formula is obtained using P(r) = r, but 
the calculations can be analytically managed for both for P(r) = —1/r and r 2 to give again (|3"3")l . The only technical 
difficulty is the inversion of the function z — W(x)x a (see App. |FJ. Since negative values for the argument T of Wo 
must be in the range [— 1/e, 0], only a limited set of values for TV are allowed. Good approximations for the critical 
values of k can be obtained by imposing e n ,i = in ((33)) @. 



A. AFM with P(r) = -1/r 

For P(r) = —1/r, we find vo = iku^ with Wo(uo)uq = —N 2 / (4k) and N = n + l + 1. AFM approximate eigenstates 
are hydrogen-like states (|C3[) with 

' 2 2 W Q (T) v ; 

In this case, g(y) = —ke 1 ^ with y < 0. Since the function g"(y) = — /ce 1 / J '(l/?/ 4 + 2/y 3 ) has not a constant sign for 
y < 0, nothing can be said a priori about the variational character of the eigenvalue (|33p with N = n + 1 + 1. But, it 
is possible to compute AV(r) = V(r, v ) — V(r) = k[e" x — r\e~ Ta jx + e~ r °(r — 1)] with r = — 3W {T). A careful 
analysis shows that AV(r) is always negative if r < 1. This implies that e y , are surely lower bounds of the exact 

energies if n + I + 1 < \Jkj (2e). Mean values (Hy; n, Z|r 2 |Hy; n, I) and (Hy; n,l\p 2 \Hy; n, 1} are then computed using 
(IC6l) and (JCZl) with (HQ). 



B. AFM with P(r) = r 2 

For P(r) = r 2 , we find vo = k/(2u ) with Wo(u )u^ 1/4 = (2/V 2 /fc) 1 / 4 and N = 2n + I + 3/2. AFM approximate 
eigenstates are hydrogen oscillator states (|D3[) with 



In this case, = —ker^ with y > 0. The function (?"(y) = —ke~^v(l/y+l/y 3 / 2 )/4 being negative in its domain, g 
is concave and e^i are upper bounds of the exact energies. Mean values (HO; n, Z|r 2 |HO; n, I) and (HO; n, l\p 2 \HO; n, 1} 
are then computed using (ID6[) and (|D7|) with (f55|) . 
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C. Results 



Results concerning some observables derived from Hamiltonian (j3"2")l are gathered in Table HXl With P(r) — r 2 , 
formula (|33[) can give an imaginary or a positive value for the energy. When a negative energy is produced, it is an 
upper bound as expected but the agreement is poor. 

Better results are obtained for P(r) = — l/r. A negative value for the energy is produced for each allowed state, 
even if the agreement is sometimes poor. This happens for highly excited states with a near zero energy. A lower 
bound is alway obtained, though the variational character can only be guaranteed for the ground state in the cases 
k = 10 and k = 20. When the value of the energy is good, observables (r 2 ) and (p 2 ) are also quite good, as well as 
the overlap. Even if the results are not so accurate as in the previous cases, analytical informations about energies, 
observables, and the number of bound states can be obtained. 

We have studied also other observables, as for the previous cases, but the corresponding results leading essentially 
to the same conclusions, they are not presented here. It is worth noting that values of energies can be greatly improved 
(see Ref. |6[) but the variational character cannot be guaranteed anymore. 



TABLE IX: Some AFM results for the Hamiltonian (J32J): R(E) = e^j/EW,!, R(r 2 ) = (X; n,l\r 2 \X; n, l)/{n, l\r 2 \n,l), R{p 2 ) = 
(X; n, l\p 2 |X; n, l)/{n, l\p 2 \n, l), and overlap |(n, Z|X; n, I) | 2 . For each value of fe, all allowed set (l,n) are presented by increasing 
values of -E n ,;, and, for each set, results are given for P(r) = r 2 (X=HO) in the first line and for P(r) — —l/r (X=Hy) in the 



second line. 


The 


symbol " 


-" indicates that the 


energy found is 


an imaginary or 


a positive numb 


ier. 




k 


I 


n 




X 


R{E) 


R(r 2 ) 


R{p 2 ) 


<n,Z|X;n,Z)| 2 


5 








-0.550 


HO 


0.245 


0.807 


1.072 


0.969 










Hy 


1.531 


0.907 


1.271 


0.993 


10 








-2.182 


HO 


0.673 


0.825 


1.123 


0.979 










Hy 


1.350 


1.136 


1.087 


0.989 




1 





-0.334 


HO 


















Hy 


1.370 


0.738 


1.346 


0.973 







1 


-0.070 


HO 


















Hy 


6.573 


0.263 


3.336 


0.465 


20 








-6.624 


HO 


0.826 


0.852 


1.117 


0.985 










Hy 


1.243 


1.275 


0.996 


0.979 




1 





-2.715 


HO 


0.647 


0.876 


1.068 


0.975 










Hy 


1.242 


1.005 


1.117 


0.992 







1 


-1.426 


HO 


















Hy 


2.365 


0.631 


1.450 


0.666 




2 





-0.431 


HO 


















Hy 


1.251 


0.773 


1.271 


0.974 




1 


1 


-0.163 


HO 


















Hy 


3.304 


0.429 


2.347 


0.625 







2 


-0.009 


HO 


















Hy 


62.1 


0.066 


11.90 


0.043 



VI. CONCLUDING REMARKS 

The auxiliary field method, which is strongly connected with the envelope theory pd| - [l6j . is a powerful tool to 
compute approximate closed-form energy formulas for eigenequations in quantum mechanics The idea is to 
replace the potential studied V(r) by a linear transformation of a well known potential P(r), the coefficients of the 
transformation being state dependent. This method was already successfully applied for a great variety of problems 
5-9]. It can also provide approximate analytical eigenstates of the Hamiltonian considered when a closed-form for 
the energy formula can be found. From a practical point of view, eigenstates for a central Hamiltonian can be 
approximated only by hydrogen- like (P(r) = —l/r) or harmonic oscillator (P(r) = r 2 ) wavefunctions which are 
properly scaled in order to match at best the studied eigenstates. 

In this paper, the quality of this approximation is tested for a Schrodinger equation with various potentials. For the 
linear potential, we found that very good results can be obtained for the wavefunction and a set of observables for the 
choice P(r) — r 2 . The typical error on energy and other mean values is around 5-10%. The error is quite independent 
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of the quantum numbers and stay constant for large values of n. Quite good results can also be obtained for the 
logarithmic potential, which is a confining interaction as the linear one. With the choice P(r) = — 1/r, analytical 
results have been found for the exponential potential but with an accuracy which is not as good as for the two previous 
cases. The main problem with this interaction is the existence of a finite number of bound states whose the most 
excited ones can have a near zero energy. 

An eigenvalue equation can also be solved within the variational method by expanding trial states in terms of 
special basis states. The correct asymptotic tail can be well reproduced if the basis states are well chosen. With this 
method, a matrix representation of the Hamiltonian is obtained and the solutions are computed by diagonalizing this 
matrix: M upper bounds of the energies are determined with the corresponding M states, where M is the order of 
the matrix. A very good accuracy is possible if M is large enough. Even for M = 1, the accuracy can be better 
than the one provided by the AFM [9]. However, if one is interested in closed-form results, the matrix elements must 
have an analytical expression and the number of computed states M must be limited to 4. But even for M = 2, the 
eigenvalues can have very complicated expression, not usable in practice. So, the variational method can only provide, 
at best, a very limited number of eigenvalues and eigenstates with an analytical form. 

The WKB method is also a popular method to solve eigenvalue equations [l|, HH- In principle, it is only valid for 
high values of the radial quantum number n, but it can sometimes yield very good results for low- lying states p4| . 
Wavefunctions are not necessary to compute some observables [24j but they can be determined for arbitrary value of n 
with this method. An advantage is that their asymptotic behavior can be correct but, unfortunately, the WKB method 
is mainly manageable for S-states. Indeed, for I ^ 0, the interaction V(r) must be supplemented by the centrifugal 
potential, which complicates greatly the integrals to compute. Moreover, these wavefunctions are piecewise-defined 
whose different parts must be connected properly at the turning points. So they are not very practical to use. 

Within the auxiliary field method, if the problem studied is analytically manageable, a state can be determined with 
the same calculatory effort for any set of quantum numbers (n, I). So this method is very useful if one is interested in 
obtaining analytical information about the whole spectra, wavefunctions, and observables of a Hamiltonian without 
necessarily searching a very high accuracy. Moreover, a AFM state can be used as a trial state to recompute the 
energy, like in the variational method. But, in this case, an upper bound can be guaranteed only for the ground state. 

It is shown in this work that the selection of the potential P(r) is crucial to obtain good results. For a linear 
potential, P(r) = r 2 is the best choice. But P(r) = —1/r provides much better results for an exponential potential, 
while the choice of P(r) is not so crucial for the logarithmic interaction. Other kind of potentials could be considered 
within the AFM. It could also be interesting to look at Hamiltonians with relativistic kinematics. 



C. Semay thanks F. Buisseret and F. Brau for useful discussions. He also thanks the F.R.S.-FNRS for financial 
support. 



In this appendix, we recall general properties for observables that are studied in this work. A given state is denoted 
\E) or |n, I) where n is the radial quantum number and I the orbital angular momentum. In order to simplify the 
notations, the mean value in, l\Q\n, I) of an observable Q is sometimes simply denoted (Q), and the abbreviation 
\n) = \n,0) is also used. 

Overlap with the ground state 

Let H be a Hamiltonian and \Ei) one of its eigenstate with energy Ei. The Eckart bound Be [25[ gives a lower 
bound on the overlap between the state \Eq) and a trial state \cf>). This bound necessitates the knowledge of the exact 
values of Eq and E\ . But it is possible to replace Be by a new bound B' E less accurate by using upper Ef and lower 
E[ bounds of a given level Ei. This gives 
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Appendix A: Observables 
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Mean values of even powers of the momentum 
Let us consider a Hamiltonian of type 



(A2) 
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where V(r) is a central potential with r — \r\. Using the Virial theorem [2(| and the hermiticity of H in the 
computation of (H 2 ), one finds 

(p 2 ) =2m(E-(V)) =m{rV), (A3) 

(p 4 ) = 4m 2 (E 2 - 2E(V) + (V 2 )) , (A4) 

with the notation Q' = dQ/dr and E = (H). Simplifications occur if the potential is a homogeneous function of 
degree A (V(f3r) = (3 x V(r)) 7 in which case r V — XV, so that these equations reduce to 

<rt-2±E and iy ) = J L- 2 E, ,A5) 

<J>*} = I" 2 (^f E * + < l ' 2 >) ■ (A6) 

This is in particular the situation for a power-law potential. 
Wavefunction at the origin 

The square of the modulus of the wavefunction ip n ,i('r) = (r\n,l) at the origin for a I = state (for I ^ 0, it 
vanishes) can be computed with a formula given in Ref. [20( : 



Generalized virial theorem 

If A is an arbitrary operator, it follows that 



IVyo(o)| 2 = -(y')- (A?) 



([H,A})=0, (A8) 



due to the hermiticity of H. If p r is the radial momentum (with [r,p r ] = i and h = 1) and f(r) an arbitrary function, 
the computation of (|A8[> with A = p r f{r) yields to the following relation, called the generalized virial theorem, 

(2J3/'(r) - F'(r)/(r) - 2F(r)/'(r) + ^ - (—)') = 0, (A9) 

\ 4m mr \ r J 

One recovers the usual virial theorem (see (IA3p ~) for the special choice /(r) = r. But new interesting relations can be 
obtained for other choices. In particular, if f(r) = r s+1 the previous equation becomes 

2(8 + l)E(r s ) - 2(s + l)(r s V(r)) - (r s+1 V'(r)) + — (s 2 - 1 - 41(1 + 1)) (r s - 2 ) = 0. (A10) 

4m v ' 

In addition, if V(r) = sgn(A)ar A for A ^ 0, (|A10[) reduces to 

2{s + l)E(r s ) - sgn(A)a(2s + A + 2)(r x+s ) + -i- (s 2 - 1 - 41(1 + 1)) (r s - 2 ) = 0. (All) 

4m 

This recurrence relation is particularly useful to compute (r k ) mean values for potentials with integer power, choosing 
s as an integer. 



Appendix B: Exact solutions for the linear potential 

The eigensolutions with I = of the Hamiltonian 

p 2 

H=— + ar (Bl) 
2m 

are analytically known in terms of the Airy functions Ai [26j |. The eigenenergies can be written in terms of the 
(negative) zeros a„ of this function, namely (h = 1) 

( a 2 \ 1/3 

E = - - a„, (B2) 
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and the normalized eigenvectors ip n o{ r ) = ('"l^) are given by 

V>„o(r) = J-^l"^ Ai f(2ma) 1/3 r + a n ). (B3) 
\/47r Ai r V / 



/47T |Ai'(a„)| r 

An approximate form for a n is given by [26j 

«» = -/3n (l + J/?" 3 - + O^ 9 ) ) will, ><„ 



3tt 

— [n+- 
2 V 4 



2/3 



(B4) 



the series converging very rapidly with n. At the origin, the square modulus of the wavefunction reduces to 



hAn,o(0)| 2 = ^. (B5) 

Remarkably this quantity does not depend on the radial quantum number. This property is specific to the linear 
potential. 

Mean values (r k ) can be computed by performing directly the integrals or by using (|A11|) . One obtains 

/ I I \ 2 j (Xn. I 

(n\r\n) 



3(2ma) 1 /3' 

8U.J2 

(n\r \n) 



(n\r 3 \n) 



15(2ma) 2 / 3 ' 
16|a n | 3 + 15 
35(2ma) 

^n| 4 + 

315(2ma) 4 / 3 



/ i 4 1 \ 8|a n | 4 + 25|a„| 

(n\r \n) = 16 . . (B6) 



Using (|A5I) and (|A6|) . it is easy to show that 



(n\p 2 \n) = (2H 2/3 ^, 

12 



H;p 4 |7i) = (2ma) 4 / 3 ^. (B7) 



Appendix C: Hydrogen-like system 

The eigenenergies of a hydrogen-like Hamiltonian 0, [l^J 



are given by {K = 1) 



and the normalized eigenvectors by 



H = f-~, (CI) 
2m r 



E = -T7^T- with JVjry=n + l + l, (C2) 



(C3) 



with 7„; = 77/(71. + Z + 1) and 77 = toza is a Laguerre polynomial and a spherical harmonic. At the origin, the 
S-states are such that 



Mean values (r k ) can be computed by performing directly the integrals or by using (|A11[) . One obtains 
k = (n + l + l) k - 1 (n + 21 + 1)! ^ , 1 ^ qrV r« (p + Q + k + 2l + 2)\ 



J2 {-iy +q cici- 



2(2 v ) k n\ ^ ~ J ~n-n(p + 2l + l)l(q + 2l+l)V 

where is the usual binomial coefficient. With L = l(l + 1), one can write also 

1\ ry 



r/ ' *V 
1 \ 2rj 

72 



" ' 8ry 4 
Using (|A5|) and (|A6|) . it is easy to show that 



(2i + l)^V 

(r) = ^( SN Hy- L )> 
N 2 

N 2 

(r 3 ) = -^-(35Nt ry + 5N 2 Hy (5-6L) + 3L(L-2)), 
1 - Hy (63AT^ + 35iV| 3/ (3-2L)+5£(3 J L-10) + 12) 



<P 2 > = * 



N 2 ' 



i . 4 8n + 2^ + 5 



(P 4 ) = V 



(2/ + l)iV4 



Appendix D: Harmonic oscillator 



The eigenenergies of the harmonic oscillator Hamiltonian [lj, |2 



are given by (h = 1) 



H = E^ + vr 2 

2m 



E = J—N H o with 2Vjj =2n + i + ~, 
V m 2 



and the normalized eigenvectors by 



with A = (2m^) 1 ^ 4 . At the origin, the S-states are such that 

l^(p)p.^5£2+M. 

7r z n! 

Mean values (r k ) can be computed by performing directly the integrals or by using (|A11[) . One obtains 
/Jts _ 1 r(n + I + 3/2) £ { _ 1)P+qcPcq T(l+p + q+(k + 3)/2) 



A* n! p ^ ' n n T(p + l + 3/2)T(q + 1 + 3/2Y 



17 



One can also write 

14I> + 3/2) 



(n\r\n) 
(n,l\r 2 \n,l) 
(n\r 3 \n) 



A mil 
N H Q 
X 2 ' 

1 8(4rt + 3)r(rt + 3/2) 
A 3 37m! 



(n, l\An, I) = Jjj Un% - 2L + . (D6) 

Using the Fourier transform of a harmonic oscillator, it is easy to show that 

(p k ) = X 2k (r k ). (D7) 

Appendix E: Overlap of eigenfunctions 

The scalar product of two radial functions R n .i(r) and R n ',i{T) of a set of orthonormal states is simply given by 
firm' ■ When one of these functions is scaled by a positive factor a, the overlap 

/>oo 

F n ,n',l(a) = a 3/2 / R n: i{x)R nl .i{ax)x 2 dx (El) 
Jo 

satisfies the following properties: 

lim F n n , i{a) = 5 nn >, 

a— >1 

\F n , n >,i(a)\ < 1. 
F n>n > t i(l/a) = F n > >ni i(a), 

lim i^,„ v (a) = lim F„. nM (a) = 0. (E2) 

a— >0 a— >oo 

The first relation stems from the definition (IE1|) . the second one from the Schwarz inequality, and the others are due 
to scaling properties. 

Using the dilation properties of the Laguerre polynomials and the various existing recurrence relations [28| . it is 
possible to compute analytically the formula -F 1 n ,n',i( a ) f° r hydrogen-like systems and harmonic oscillators. In the case 
of hydrogen-like systems, one obtains 



C'» = (N + l)\n'\ (N' + 0! (iaNN'f s ff + 7+i 



E(-D fe 



X 
k=0 



1 



4aiViV' J k\(n- k)\ (N — k + 1)\ (n' — n + k + 1)! 



x [2(N - fe)(n' -n + fc + l) + (n- k)(N -k + 1)^7 + {n - n + k){n' - n + k + l)^r) 
\ 2aJS Q{a) J 



(E3) 



with N = n + I + 1, N' = n' + I + 1, Q(a) = a./V — AT', and S(a) = aN + N'. In the case of harmonic oscillators, the 
formula is given by [291 ] 

2\ n ' ~ n 

F™ ri (a) = y/n\ n'l T(n + l + 3/2) T(n' + I + 3/2)(2a) 2n+l+3 / 2 ' 



(1 + a 2y i + n '+ l + 3 / 2 

n /i _ 2^ 

^ ' {2a) 2k k\(n-k)\(n' -n + k)\T{n-k + 1 + 3/2)' { ' 



It can be checked that both (|E4[) and (IE3|) satisfy relations (IE2[) . but the calculations are much more tedious for F H y . 
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Appendix F: Solving a special equation 

The multivalued Lambert function W(x) is defined by W^ 1 (x) — xe x which implies that W{x)e w ^ = W(xe x ) = x 
(see Ref. @ for more details). In section [V] it is necessary to solve the following equation (giving x as function of the 
parameter z): 

z = W{x)x a . (Fl) 
Using the two changes of variables x — ye v and ay = (a + this equation can be inverted to give: 

x = W~\z) = ze z if a = 0, 

x = W- 1 (In- ] = -In- if a = -1, 

\ z J z z 

x = w ~ x i^r w (^tt z1/(q+1) )) otherwise - ( F2 ) 

The domains of the functions z(x) and x(z) must be examined carefully. They depends on the branch Wo or W-\ 
considered for the Lambert function and on the value of a. 
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